Assessing structural bonding aspects of multiband 
superconductors through impurity-induced local 
lattice distortions: a case study on MgB2 



Aleksandr Pishtshev^ and Mihhail Klopov^ 

^ Institute of Physics, University of Tartu, Riia 142, 51014 Tartu, Estonia 
E-mail: aleksandr . pishtshevOut . ee 

^ Department of Physics, Facuhy of Science, Tahinn University of Technology, 
Ehitajate 5, 19086 Tahinn, Estonia 

E-mail: mihhail.klopov@ttu.ee 

Abstract. We report the results of an ab initio modeling of substitutional impurities 
such as zinc, copper and zirconium ions incorporated into the magnesium sublattice of 
MgB2. These species are of particular interest because they demonstrate the lowest- 
observed-droop of the superconducting temperature among most of the impurities 
of the same formal valences. The goal of computational studies was to gain an 
atomistic understanding of how the structural and bonding properties of the local ionic 
environment are affected by a change in the host cation. The simulations performed 
for the given set of substituents indicate that at a low doping level, the induced 
lattice distortions and additional forces are noticeable small. At the same time, the 
electron redistributions around the impurities have been found to differ drastically 
both from that around the host cation and from each other. Results of the first- 
principles calculations were used to compare and discuss, in the context of the cation 
properties, specific changes in local charge structures and in chemical bonding caused 
by the presence of the impurity. Certain trends in substituent effects were found. 
It was determined how the substitutional ion modifies the initial charge distribution 
around the host cation position and therefore, affects both the overall picture of the 
electronic states involved in the chemical bonding and the intrinsic charge-transfer 
channels. Our studies demonstrate that the robust structural stability of the boron 
sublattice of MgB2, which preserves the nearest-neighbor environment around the B 
atoms practically unchanged from an impurity incorporation into the magnesium site, 
has important functional significance ~ it supports a much weaker dependence of the 
superconducting temperature on impurity content. The present work gives insight into 
structural and electronic features of the impurity caused-processes which contribute to 
unique material properties of various MgB2-solid solutions. 
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1. Introduction 

In a theoretical analysis of impurity signatures, principally important is the fact that the 
equilibrium between the (nearest neighbors) lattice atoms, and the implanted impurity 
atom represents the structural balance between local geometry and chemical bonding, 
which depends strongly on the interatomic distances. In order to explore the dynamical 
nature of substituent effects on small length scales one can employ an approach based 
on a detailed atomistic view: it can be straightforwardly implemented within ab-initio 
numerical simulations, which are able to calculate the relative ionic displacements and 
the charge density redistribution accompanying the effect of the substitution. The 
results of calculations will then allow us both to capture the final binding geometry 
around the impurity site and to answer the fundamental question as to what extent the 
impurity atom is capable of adapting itself to the local environment and the existing 
chemical bonds. 

In this paper, we apply such integrative approach to investigate the effects 
of substitutions for several 3d transition- metal impurities such as divalent Cu(II), 
Zn(II) and tetravalent Zr(IV) in magnesium diboride MgB2. We consider MgB2 
containing substitutions of Mg as a doped model system for which we perform a 
series of the relevant numerical simulations. There are two reasons that make this 
system attractive for detailed study. First, while different properties of magnesium 
diboride doped by a variety of impurity atoms have been extensively studied in the 
past decade, the investigation of impurities in this material is still a topical subject, 
raising a number of actual issues that are worthy of special attention (for details 
see P 121 121 H El il m El m Uni HI [12], and the references cited therein). The other 
reason is that knowledge on how the arrangement of a single impurity ion influences 
on a local structural order may have useful predictive value in the rational design of 
magnesium diboride solid solutions. One of those issues, which is especially important 
to all aspects of local geometry and stability, is the influence of substitutions on a 
principal atomic configuration and the bonding properties of the cation environment. 
However, in many cases, information on the bonding situation around a substitutional 
impurity turns out to be limited, and there remain several questions about binding of 
metal elements incorporated into the place of the host Mg(II) cation. Understanding 
the latter is very essential because the response of the system on doping is not clear 
a priori. For example, it was experimentally established that the suppression of the 
superconducting properties proceeds much slower in the Zn-, Cu- and Zr-doped systems, 
as compared to compounds that contain other dopant ions of the same formal valences 
(details can be found in [H IH El |9l HH [13]). In table 1, using the data available from 
the literature, we have summarized some examples related to this observation. 

From the rate of decrease of the superconducting temperature Tc with chemical 
substitution of the original Mg with metal elements in MgB2 one can expect that not 
all the impurity properties are matched well with those of Mg. Thus, the given example 
of table 1 shows that conceptually Zn-, Cu- and Zr-containing solutions stand among 
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Table 1. Experimental lowering of the Tc value (ATc) in the MxMgi_xB2 compositions 
at substitution levels (x) in the range 0.03 < a; < 0.10. 



M Zn 


Cu 


Zr 


Ca 


Co 


Fe 


Ni 


Pb 


X : ATc 0.03 : -0.2^ 


0.03 : 0.5- 


0.05 : 0.5- 


0.03 : 0.5s 


0.03 : 2.8*^ 


0.03 : 5.r 


0.03 : 2.9' 


0.03 : 1.5j 


0.05 : 0.5'^ 


0.05 : 0.2'^ 


0.05 : 0.7^ 


0.05 : 1.2s 


0.03 : 2.0' 


0.05 : 8.5' 


0.05 : 4.8' 


0.05 : 1.8j 


0.10 : 0.2'' 




0.10 : 1.1^ 


0.07 : 1.8S 


0.05 : 3.3' 






0.07 : 2.2j 


-[a; b[Il]; -[IS]; d[n]; 




=[in];''[2I];^[II];j[ll] 











other doped systems as a particular example. This feature prompted us to examine in 
more detail the role of these substitutional ions and ion-specific effects in the context of 
impurity-induced changes in local structure, dynamics and chemical bonding. Therefore, 
the aim of the present paper is by employing ab initio calculations to study variations of 
certain lattice and charge characteristics, such as lattice distortions, interatomic forces, 
charge densities and bonding features, caused by compositional variations in MgB2. By 
comparative analysis of features of the impurity atom with those related to the host one 
we pay attention in our work to what is happening within an ionic environment in the 
elementary cell when an intensity of interionic and electronic correlations is affected by 
the replacements of a magnesium cation. 

2. Computational details and impurity simulation 

According to the literature data (e.g., [H HJ [23]), a crystalline structure of the bulk 
MgB2 consists of alternating honeycomb layers of B interleaved with closed packed 
layers of Mg. Each Mg atom is located at the center of a hexagon formed by boron; 
by donating two valence electrons it becomes a strongly charged ion with a charge 
very close to +2. Consensus on the nature of the superconductivity in MgB2 is 
based on a generic multiband model which employs the intra- and interband pair- 
scattering mechanisms. The presence of the interband channel connecting the a and 
the 71 electronic bands [211 [25], [26] gives rise to a unique property of MgB2 - two- 
gap superconductivity, which is characterized by two gaps of different sizes that are 
simultaneously closed at the same temperature Tc of the superconducting transition. 

The first-principles modeling of the effects of Mg(II) ion substitutions in the 
MxMgi_xB2 compositions, where M stands for a transition metal cation such as Zn, 
Cu or Zr, was performed using the Vienna Ab-initio Simulation Package (VASP) [271 
[281 [29| [30] in conjunction with a plane-wave-basis code and the projector- augmented 
wave (PAW) method [311 [3^ . Exchange-correlation effects were described through the 
Perdew-Burke-Ernzerhof (PBE) framework [33]. For geometry optimization, one of 
the host atoms was replaced by one M element in a 81-atoms supercell (3x3x3). 
Such arrangement shown in figure 1 corresponds to a; = 0.037 in terms of the formal 
stoichiometry of the doped system. The lattice distortion caused by the substitution 
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has the following features: (1) it is totally symmetric on the cation site, and (2) it 
involves displacements of 12 nearest-neighbor boron atoms from regular positions. For 
a central boron atom, these displacements lead in turn to local symmetry breaking. 
Full relaxation of atomic positions was made with periodic boundary conditions and 
theoretical (structurally relaxed) bulk lattice parameters. The convergence criteria for 
total energy and residual forces calculations were set to be within 1 x 10~^ eV/atom and 
1 X 10^^ eV/A, respectively. A cut-off energy of 700 eV and a Monkhorst-Pack A;-point 
grid [34j of 5 x 5 x 5 has been employed. This computational scheme was successfully 
tested by a comparison of the calculated values of the lattice constants, acaic = 3.069 
A; Ccaic = 3.504 A; Ccaic/dcaic = 1.142, with the experimental data [23], a = 3.085 A; 
c = 3.523 A; c/a = 1.142 (the deviations are 0.52% and 0.53%, respectively). 

Changes in chemical bonding upon incorporation of a dopant atom into a Mg-site 
was studied in terms of charge-density distributions. In order to characterize variations 
of the partial charges at the magnesium position upon substitution and to assign a 
specific meaning of ionic charges in the analysis, values of the Bader effective charges [3S] 
have been evaluated. Their calculations were carried out by using software [36], which 
utilizes a grid-based Bader analysis algorithm [37]. The contribution of local many-body 
polarization effects was analyzed in terms of the Born effective charges, values of which 
have been obtained by using the special extended method adopted at VASP. 

3. Results 

In tables 2 to 5, we present the main results which characterize the influence of the Zn, 
Cu and Zr substitutional impurities on displacements of the boron atoms, the electronic 
charge redistribution around a cation, and the local properties of chemical bonding in 
MgB2. These results and several consequences of magnesium replacement are discussed 
below. The main conclusions are summarized at the end of the paper. 

3.1. Characterization of local distortions 

We pay attention to the data of table 2 as they reveal certain trends. First, no significant 
changes are seen - local lattice distortions, and the induced changes in transverse ionic 
strengths caused by cationic substitution appear to be notably small. Particularly, this 
observation relates to cases with Cu and Zn. This implies that the relevant bond lengths 
are preserved, and the boron layer remains planar. Further, it is also seen from table 2 
that longitudinal (in-plane) stresses are more strongly affected by the impurities, though 
the induced changes remain relatively small. 

The overall comparison of the data of table 1 with the simulation results of table 2 
reveals a certain correlation between the variation in the magnitude of the Tc suppression 
rate (ATc/Ax) and the doping induced changes in the boron sublattice. This can be 
emphasized by observing the following difference: While in the majority of the doped 
solutions, the presence of several percents of divalent impurity metal atoms provides 



5 



Table 2. Calculated lattice distortions and changes of interatomic forces in the boron 
layer caused by substitutional impurities in the magnesium diboride'*. 



Host/Imp. 
cation 


B out 
(A) 


B in 

(A) 


force components 
along c-axis in a-b plane 
(eV/A) (eV/A) 


Mg 














Zn 


+0.0045 


-0.0054 


+0.020 


-0.089 


Cu 


-0.0002 


-0.0128 


-0.014 


-0.234 


Zr 


-0.0063 


+0.0159 


-0.013 


+0.260 



^Denoted: 'B out/B in' - boron out-of/in-plane displacement as shown in figure 1; the 
minus sign indicates inward directions. 

already a considerable suppression of Tc, there are substitutions such as Zn, Cu, or even 
tetravalent Zr low concentrations of which slightly influence both the superconducting 
order and the specific positions of the boron atoms. The reason for this discrepancy 
seems to be related primarily to the strong 2D B-B linkage dominating in the boron 
layer [38] which, as we have shown above, tends to keep the existing structural stability 
upon a change in the host cation. Considering Mg as a reference cation, another 
interesting result is that one may say according to [3l [39] that the incorporation of 
just Zn, Cu or Zr into the Mg-site has little effect on actual interband pairing channels. 
Correspondingly, in these cases, the protection function of the covalent boron linkage 
cannot be markedly influenced (at low doping levels) what may provide a much weaker 
dependence of the superconducting temperature Tc on impurity content. 

3.2. Characterization of dopant local bonding 

In order to get further insight into why the substitutional Zn, Cu and Zr cause only 
minor changes in a planar geometry of the covalent boron layer, we undertook an 
investigation of how local properties of chemical bonding are affected by the overall 
electronic structure and elemental characteristics of each of the considered impurity 
ions. For characterization of impurity-induced changes in chemical bonding via the 
differences of the electron states between the host and foreign cations, we performed 
an extensive comparative analysis of the relevant charge partitioning schemes. Figure 2 
shows charge density maps calculated for MgB2 and all three doped cases; the maps are 
projected onto the (100) plane and contain the central Mg or M = Zn, Cu, Zr cations, 
and surrounding B atoms. Figure 3 is instructive in another way: for each cation 
considered, it shows the line charge density distributions along the Mg-B and M-B 
directions, respectively. Computed values of the Bader effective charges (Qb) presented 
in table 3 reflect the topology features of the underlying electron-density distribution 
emerging as a direct consequence of the relevant arrangements. Comparison of figures 
2, 3 and the data of table 3 for the reference case of the Mg(II) ion reproduces the 
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Table 3. Some charge characteristics related to the host cation and its substituents 
in the magnesium diboride. The Pauling electronegativities are taken from [3D]. The 
quantities Qb(M) denote Bader effective charges calculated from electronic densities 
(in units of |e|). 





Mg 


Zn 


Cu 


Zr 


Electronegativity difference 
values with respect to boron 


0.73 


0.39 


0.14 


0.71 


Qb(M) 


+2.0 


+0.69 


+0.42 


+2.68^ 



comparison shows that this value appears to be equal to that of +2.67 obtained 
from the self-consistent LAPW calculations for the bulk zirconium diboride ZrB2 [41]. 

well-established result [21 [23lll2l [l3l HHIIH] that Mg is an ionically-bonded cation (with 
an admix of some weak covalency) in a stable divalent state. 

Analysis of cases where Zn or Cu substitutes for Mg revealed three significant 
differences. First, as depicted in table 3, the Bader effective charges are quite different 
from the formal chemical charge +2. Two other differences are very clearly delineated on 
figures 2 and 3; common to both ions they are: a depletion of the electron charge along 
the metal-boron bonding direction, and a large degree of electron charge accumulation 
of fairly good spherical form around the nuclei. Correspondingly, the physical picture of 
the outer electronic structure of Zn and Cu incorporated into magnesium diboride can 
be drawn as follows: (i) the 3d states appear to be completely filled (what provides 
increased electron density at the substituent, causes the observed spherical charge 
distribution around it, and induces the smaller effective charge transfer, in comparison 
with magnesium), (ii) zinc and copper cations acquire fractional quantities of the 4s 
electron populations, about 1.31e and 0.58e, respectively, (iii) the partially filled 4s states 
constitute the resulting outermost orbitals, and (iv) the d-electron covalency effects are 
extremely low because the filled 3d levels of the substituent are shielded by its 4s states. 
Remarkably, smallness of the effective charges relative to the nominal chemical valence 
+2 is an indication of the fact that the substitutional cations have retracted the certain 
amount of the delocalized electrons that participate in metallic bonds. The reason why 
just the electrons shared by the boron atoms are excepted is that the significant changes 
in an equilibrium planar geometry of the covalent boron layer are absent. Furthermore, 
of note is another observation from table 3 that the local loss of charge exchange activity 
is directly connected with the drop in the electronegativity difference with respect to 
boron in the series of elements Mg, Zn, Cu. 

It may be concluded thus that both substitutions give rise to a negative tendency 
which manifests itself (i) in increasing the electron density around outer orbitals of 
the Zn or Cu cations, and (ii) in decreasing the positive charge at the cation center 
up to the values of QB(Zn) = +0.69 and Qb(Cu) = +0.42, respectively. This leads 
to partial neutralization of the overall charge on the substituents, makes the initial 
charge transfer channels much weaker and therefore causes the significant reduction of 
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the polar contribution into the cationic nature of the substitutional ion. As a result, 
the initial, predominantly ionic character of the local environment is degraded. The 
variances of the Bader effective charges from the nominal +2 may therefore be regarded 
as a degree of Zn or Cu "underbonding" that should have a certain effect on physical 
properties of MgB2 doped by these elements. For example, since the interband channel 
is unaltered (because of minor changes in a planar geometry of the covalent boron 
layer), the experimentally observed smallness of suppression of the superconducting 
temperature is perfectly understandable via influence of the above-mentioned retraction 
of "previously expelled electrons" on the density of electronic states at the Fermi level. 
This result correlates well with the conclusion made in |16] for a series of Al-Li codoped 
MgB2 systems that the superconductivity is much more affected by the relevant lattice 
distortion induced by impurity ion substitution than by the band filling. 

The trend towards electron density accumulation and a weakening of charge transfer 
channels breaks down for the case of the substitutional Zr in the structure of MgB2. 
Although Zr possesses almost the same electronegativity difference as Mg , the calculated 
value of the Bader effective charge, QB(Zr) = +2.68, is considerably different from 
the formal valence 4+. From chemical point of view, this implies that the zirconium 
cation in the MgB2 lattice tends towards a lower oxidation stat^. In contrast to zinc 
and copper ions, Zr attempts to expand definite amounts of charge density from itself 
toward an interstitial area. Analysis of figures 2 and 3 strongly suggests the moderate 
enhancement of the positive charge on the cation which takes place due to delocalization 
of its 4d electrons caused by partial donation of 0.68 additional electrons. Such a 
delivery of a partial charge (along the metal-non-metal charge-transfer configuration) is 
well correlated with a chemical readiness of the zirconium cation to give up the valence 
4d electrons in order to move back to the predominant oxidation state of 4+ and to 
activate the covalent bonding with neighboring atoms. Consequently, as compared with 
the parent compound MgB2, an ionization increase Mg^+ — > Zr^-^^+ on a substitutional 
cation site and the presence of appreciable covalency in the same region may give rise 
to a tendency towards dielectrization of the system. This observation is confirmed by 
first-principles studies of the bulk compound ZrB2 [ISl SH SS] , which show the presence 
the strong hybridization between the Zr 4d and the B 2p states as well as drastic fall of 
the density of electronic states at the Fermi level. 

Thus, we can conclude that electronic charge distribution at the Mg position in 
the bulk magnesium diboride is very sensitive to the cation replacement. Firstly, the 
substitution for an impurity ion of the same nominal charge +2 leads to a substantial 
charge redistribution around the cation. Secondly, the change in the chemical bonding 
character depends on the type of metal ion and its electronegativity versus boron, so 
that among substitutional ions considered two different trends can be distinguished. All 
this underlines in the light of the bonding situation that ion substitution not only always 
affects the existing chemical bonds but also that the resulting charge rearrangement is 

J This point suggests the presence of some mechanism of valence reduction upon substitution, as can 
be compared with 4d^ configurations of low-valent Zr(II) in complex sheet and cluster structures [47] . 
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exhibited quite differently for the cations considered in the present work. Moreover, 
in the context of superconductivity, these features can seriously limit the targeting 
efficiency of a number of potential substituents or even make them unfavorable towards 
the enhancement of superconducting properties. 

3.3. Pressure effects 

To complete the results obtained for impurity-induced changes in the cation charge 
states, it should be interesting to study how these changes proceed under pressure. For 
this purpose, the Zr- and Zn-doped case were considered. We employed the same set 
of calculations as used throughout this work and fully relaxed a doped system. The 
resulting pictures of the charge density maps are shown in figure 4. The comparison 
of the charge distributions with those shown in figure 2 reveals no significant changes 
between them. The most remarkable difference is a slight increase of the Bader effective 
charge of the Zr cation - from +2.68 at zero pressure to +2.84 under pressure of 
16.5 GPa (what evidently indicates the above-mentioned oxidation tendency to the 
tetravalent state). This implies that the charge state of Zr remains relatively stable 
with compression until at least 16.5 GPa. 

The Zn substitutional cation also demonstrates the similar tendency to maintain its 
charge state under pressure. For example, the calculated change in the Bader effective 
charge under compression of 14.7 GPa, AQB(Zn) = +0.04, shows that due to contraction 
of the filled 3d levels, valence electrons do not favor extra oxidation. As a result, the 
overall charge configuration of Zn in magnesium diboride is essentially unaffected by 
pressure. 

3.4- The Born effective charges 

Another possible microscopic reason of lattice distortions stems from the fact well-known 
in the theory of local phase transitions that distortion of local geometry around an 
impurity may be caused by sufficient strength of the purely dynamical part of the charge 
transfer (e.g., [50| IST]). Generally speaking, significance of dynamical charge transfer 
processes generated by vibronic mixing of the proper electronic states is one of most 
important characteristic features of polar crystals [521 |53l [Ml |55] . On the atomic level, 
this feature is associated with many-body polarization effects that contribute to ionic 
dipole polarization, and it reflects a mixed ionic(polar)-covalent character of the relevant 
chemical bonds [55]. Since the latter can be interpreted in terms of the Born dynamic 
effective charges with respect to some reference (nominal) ionic value (e.g., [561 155]). 
the calculations of the principal values of the Born dynamic effective charge tensor Z* 
can be employed as a test of how electronic effects of the many-body nature translate 
the polar contribution into ionic chemical bonding. It is also of note that although the 
idea of calculation of the Born dynamic effective charges for a conducting material may 
seem senseless at first glance, it can be, due to logic behind the ionic bonds and charge 
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Table 4. Principal values of the Born effective charge of Mg(II) in magnesium diboride 
(in units of |e|). The diagonal form with Z*x(Mg) — Z*y(Mg) reflects hexagonal 
symmetry of the tensor Z*(Mg). The quantity Z*(Mg) represents the average over 
crystal axes: Z*(Mg) = (l/3)(Z*,(Mg) + Z* (Mg) + Z:,(Mg)). 



Z*x(Mg) 


z;,(Mg) 


Zz*z(Mg) 


Z*(Mg) 


+ 1.87 


+1.87 


+2.31 


+2.02 



transfer, reasonably useful for a mixed bonded solid, indicating how specific pieces of a 
composite picture of chemical bonding may be distinguished. 

In table 4, we have presented the converged values of the Born effective charge 
for the magnesium cation in MgB2 (due to metallic properties of MgB2, we focused 
on dynamic charges only for cations and not for boron atoms). It may be assumed 
that the numerical convergence to stable values of table 4 is guaranteed by the nearly 
ionic character of the Mg sublattice which, in the spirit of layered superconductivity 
approach [571 158] . can be figuratively imagined as a stack of the intrinsic insulating 
layers. Calculations have also predicted principal values of the electronic dielectric 
constant for the bulk magnesium diboride: = = 22.45, = 5.77. The theoretical 
value eoo = 16.9, evaluated as an average eoo = (l/3)SieJ^, is quite comparable with the 
experimental estimate of 11.9 [8] thus proving the consistency of the present analysis. 

The tensorial nature of the Born effective charge provides subtle information 
concerning differences in directions and strengths of dynamical processes related to 
the ionic sites and caused by the polarization effects. As seen by comparing the values 
of table 4, there exists a certain anisotropy, as is indicated by the observation that in 
the directions perpendicular (the a6-plane) and parallel to the c-axis the components 
of the dynamic charge differ markedly from the formal value of 2+. This suggests 
that, although on the average Mg remains quite well ionized in the full A^-particle 
ensemble (Z*(Mg) = +2.02), (i) the attached effective potential, which gives rise to 
polarization force fields acting in a local environment around magnesium, appears non- 
spherical, and (ii) due to deviation from the nominal ionic charge some contribution of 
the pure dynamic covalent nature in the magnesium-boron bonding should be assumed. 
The latter is entirely consistent with previous studies of dynamical covalency effects in 
MgBa |59]. 

In order to analyze how the change in the host cation influences the attached 
effective potential, we carried out the relevant calculations of the Born effective charges 
for the substitutional impurities we are considering here. In contrast to the previous 
case of magnesium, the computational method failed to give uniform results consistent 
with hexagonal symmetry. In the physical sense, the lack of numerical consistency 
is not discouraging, and should be considered as an additional confirmation of the 
drastic changes of a charge state which arise around the host cation position due to 
its replacement. In other words, impact of the substitutions with such unavoidable 
outcomes as falling of ionic character of chemical bonding, and the expansion of a 
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Table 5. The averaged principal values of the Born effective charge tensor for the 
foreign cation M — Zn,Cu,Zr in magnesium diboride (in units of |e|). Z*''*^^) (M) = 
(l/3)(Z*l'^^s'(M) + Z;("^*^)(M) + 





Z;r^)(M) 


z;^«)(M) 


Zzl'^^s^(M) 


Z*(avg)(-^) 


Zn 


+1.6617 


+ 1.6617 


+2.25 


+1.86 


Cu 


+1.1378 


+ 1.1345 


+1.85 


+1.38 


Zr 


+1.4018 


+ 1.4630 


+2.26 


+1.71 



metallic environment beyond the range of the boron layers [ISl SB] niake it impossible 
to get from the method of calculations the target values of dynamic charges for the 
foreign cations. To obtain some approximate estimates of the Born effective charges, we 
have constructed a smooth procedure based on averaging the obtained results over the 
supercell, so that all the cations employed in the calculations are taken into account. 
This can be done since overall charge balance is maintained. Due to electroneutrality 
of the simulation supercell, one can then expect that summation eliminates off-diagonal 
contributions, suppresses discontinuities, random perturbations and other artificial 
symmetry effects and thus will return the weighted-average components of the Born 
effective charge tensor in the original diagonal form. 

In table 5, we have presented the numerical values of the dynamic charges for 
Zn, Cu and Zr substitutional cations as averages over the supercell. It is seen that 
the remaining uncertainties on the averaged values of charges are really negligible, so 
that this type of approximation is quite reasonable. Nevertheless, one additional aspect 
which should be discussed here is the sensitivity of the retrieval accuracy within such 
a calculation procedure. First of all, we note that all off-diagonal contributions are 
exactly eliminated. On the other hand, the required identity of the diagonal terms 
(Zxl^^^^ = Zyy*^^^^) is restored at different accuracy levels: for zinc and copper, we get a 
more accurate retrieval than for zirconium. We believe that this result is not surprising 
because some discrepancy in the case of zirconium suggests that due to the expansion 
of metallic bonds into the cationic environment, the effect of charge fluctuations is not 
completely eliminated after summation over the supercell. 

The results listed in table 5 allow us to characterize the specific features of chemical 
bonds in terms of the Born effective charges as follows, (i) A dynamic component in 
resulting chemical bonding becomes smaller upon cation substitution (as compared with 
that of Mg), and (ii) this partial loss in initial dynamic covalency is caused by the 
relevant charge redistribution around the cation position. At the same time, as seen 
from table 5, a rate of an anisotropy in the dynamics of valence charge distribution 
remains relatively unchanged in the cases of zinc and zirconium. Moreover, of all three 
components just the z component of the Born effective charge, which is particularly 
interesting for dynamical covalency effects in magnesium diboride, tends to be close to 
the c-axis value of the Mg cation and therefore appears, in comparison with in-plane 
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X, y components, to be more robust to alteration of charge densities. Following the 
same consideration as in [60], this leads us to an understanding regarding why lattice 
distortions induced by Zn, Cu or Zr ions are so small: the explanation is that the ionic 
strength changes induced by variations of dynamical charge transfer along the c-axis are 
not large enough to drive the relative boron displacements. 

4. Conclusions 

In order to understand the leading factors that govern adjustments of the lattice 
positions during accommodating the foreign cation in MgB2, we investigated the induced 
deformations (lattice distortions) in the local geometry of the boron sites. The first- 
principles calculations were performed for three different substitutional cations such 
as divalent Zn(II) and Cu(II), and tetravalent Zr(IV). These species are of particular 
interest because they demonstrate the lowest-observed-droop of the superconducting 
temperature among most of the impurities of the same formal valences. Inspection of 
the relaxed lattice structures together with the results of the performed calculations 
allowed the following conclusions to be drawn: (i) At a low doping level, the Zn, Cu 
and Zr substitutions cause only negligible changes in the structural arrangement of the 
neighboring boron atoms, (ii) Such an insignificant impact on the local geometry of the 
boron layer is mostly attributed to two aspects. The first is the observation that the 
existing in-plane B-B covalent bonding is strong enough to keep effectively the boron 
atoms non-shifted under substitution of magnesium with Zn, Cu or Zr. The second is 
that dynamical covalency processes along the cation-boron bonding tend to be tolerant of 
changes in the host cation and, consequently, cannot create a channel of instability in the 
boron layer. Both aspects are important in view of the fact that the structural balance 
between local geometry and chemical bonding are related to the electron structure of 
the impurity and to the bonding connections carried out by the impurity, (iii) A role 
of small deformations of the local geometry of the boron layer induced by the impurity 
ions may be insufficient to influence alone on superconductivity of the host system, (iv) 
Generally, our results emphasize that MgB2 is the perfect material having a unique 
combination of an ideal planar layer structure, the predominantly ionic character of the 
cationic sublattice, and the high degree of chemical bonding stability. 

In order to investigate the characteristic measures of impurity activity, such as 
charge states and the impurity-induced changes in the chemical bonding, we performed 
the comparison analysis of the underlying charge-transfer processes. It has been 
discovered that the replacement of Mg(II) by divalent Zn(II), Cu(II) or tetravalent 
Zr(IV) influences substantially the character of the charge transfer, so that the 
distributions of the charge around the impurities and degree of ionicity differ drastically 
both from that around the host cation and also from each other. It has been shown that 
the different cationic nature of the substituents dictates two entirely opposite directions 
along which the induced changes in the local structure take place. First, the change in 
the host cation leads to density contraction - the substitutional Zn or Cu demonstrate 
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the considerable increase of the overall charge density due to retaining its valence 
electrons and the small amount of charge transfer to the boron layers. Particularly, 
a larger accumulation of a charge and the least charge transfer are observed for the 
copper ion. It was argued that this process retracts the certain amount of the delocalized 
electrons from metallic bonds what causes the observed decrease in the superconducting 
temperature values as the substitution level grows. Second, the change in the host 
cation leads to charge delocalization - no charge accumulation is observed for Zr. The 
Zr incorporation into the Mg site of MgB2 tends (a) to preserve the initial charge 
transfer channels, and (b) by partial donation of 4d electrons to add some amounts of 
ionicity and covalency into the local cation-anion coupling. It was suggested that such 
increasing of covalency character does not favor superconducting properties of the host 
material. The stability of charge states related to the impurity cations was tested against 
pressure: the pressure dependence of the charge transfer degree was found to be very 
weak. The theoretically observed effective charges indicate that the substitutional ions 
do not possess their formal valence (oxidation state) , and are characterized by different 
levels of electron delocalization; zirconium as an electron donor is indeed stronger than 
magnesium. As a consequence, by comparison with less positively charged zinc and 
copper, the behavior of zirconium remains much more tolerant (affine) to the usual 
ionic distribution associated with the host magnesium cation. 

Summing up, one can say that the induced local lattice distortions and changes of 
the charge transfer degree may serve as structurally sensitive prognostic parameters for 
signaling how the superconductivity in MgB2-solid solutions is affected by impurities, for 
clarifying the functional roles and efficiency of such various factors as a geometry of the 
boron layer, the electronegativity difference, ionic character of the local environment, 
relative affinities toward magnesium and boron, etc., and on this basis for facilitating 
more precise theoretical predictions in synthesis and characterization. 
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Figure Captions: 

Figure 1: A bulk 81-atoms supercell modeling the MxMgi_xB2 composition for 
X — 0.037. The MgB2 cell containing the impurity cation is sketched. In-plane and 
out-of-plane boron displacements are indicated by arrows. 

Figure 2: Charge density maps of (a) MgB2, (b) Zno.o37Mgo.963B2, (c) 
Cuo.o37Mgo.963B2 and (d) Zro.o37Mgo.963B2 projected onto the (100) plane. 

Figure 3: Comparison of the line charge density distributions calculated along the 
cation-anion direction. The dashed, dot-and-dash and solid curves represent 
Zr-substituted, Cu-substitutcd and Zn-substitutcd MgB2, respectively. The inset 
shows the corresponding profile for the non-substituted MgB2. 

Figure 4: Charge density maps of (a) Zr-substituted and (b) Zn-substituted MgB2 
projected onto the (100) plane. Calculations were made for a pressure of 16.5 GPa and 
14.7 GPa, respectively. 
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Figure 1. A bulk 81-atoms supercell modeling the MxMgi_xB2 composition for 
X — 0.037. The MgB2 cell containing the impurity cation is sketched. In-plane and 
out-of-plane boron displacements are indicated by arrows. 
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Figure 3. Comparison of the line charge density distributions calculated along the 

cation-anion direction. The dashed, dot-and-dash and solid ciirves represent Zr- 
substitutcd, Cu-substitutcd and Zn-substitutcd MgB2, respectively. The inset shows 
the corresponding profile for the non-substituted MgB2. 
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Figure 4. Charge density maps of (a) Zr-substituted and (b) Zn-substituted MgB2 
projected onto the (100) plane. Calculations were made for a pressure of 16.5 GPa and 
14.7 GPa, respectively. 



